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ABSTRACT 



o 

Cu . Aims. To model the chemical evolution of manganese relative to iron in three different stellar sys- 

tems: the solar neighbourhood, the Galactic bulge and the Sagittarius dwarf spheroidal galaxy, 
and compare our results with the recent and homogeneous observational data. 
C^ [ Methods. We adopt three chemical evolution models well able to reproduce the main properties 

^"^ ' of the solar vicinity, the galactic Bulge and the Sagittarius dwarf spheroidal. Then, we compare 

different stellar yields in order to identify the best set to match the observational data in these 



t"*"- ' systems. 

^Ti ' Results. We compute the evolution of manganese in the three systems and we find that in order to 

c7^ ■ reproduce simultaneously the [Mn/Fe] versus [Fe/H] in the Galactic bulge, the solar neighbour- 

hood and Sagittarius, the type la SN Mn yield must be metallicity-dependent. 
Conclusions. We conclude that the different histories of star formation in the three systems are 
?— ^ ' not enough to reproduce the different behaviour of the [Mn/Fe] ratio, unlike the situation for 

[ff/Fe]; rather, it is necessary to invoke metallicity-dependent type la SN Mn yields, as originally 
suggested by Mc William, Rich & Smecker-Hane in 2003. 

Key words, nuclear reactions, nucleosynthesis, abundances - Galaxy: abundances - Galaxy: 
evolution 



1. Introduction 

McWilliam, Rich & Smecker-Hane (2003) considered the abundance trend of Mn in three contrast- 
ing stellar populations: a sample of Galactic bulge K giants, stars in the solar neighbourhood and 
stars in the dwarf spheroidal galaxy Sagittarius (Sgr dSph). They found that the trend of [Mn/Fe] 
in the Galactic bulge follows more or less the relation of the solar vicinity, but most stars of Sgr 
dSph show [Mn/Fe] deficient by ~ 0.2 dex. They concluded that to explain this result one has 
to assume that the yields of Mn both from type II and type la SNe should depend on the initial 
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stellar metallicity. This conclusion was at variance with a suggestion by Gratton (1989) that Mn is 
overproduced in type la SNe, relative to type II SNe. Gratton (1989) observed that [Mn/Fe] versus 
[Fe/H] in solar vicinity stars show the opposite behaviour of the [a/Fe] ratios: in fact, [Mn/Fe] 
~-0.4 dex and constant for stars with [Fe/H]< -1.0 dex and increases up to the solar value for 
stars with [Fe/H]> -1.0 dex. As it is well known, first Tinsley (1979) and then Greggio & Renzini 
(1983) and Matteucci & Greggio (1986) explained the trend of the a-elements (e.g. O, Mg, Si, Ca, 
Ti), relative to Fe, as due to the different roles played by type la and II SNe in galactic chemical 
enrichment. In particular, the higher than solar [a/Fe] ratio in Halo stars is interpreted as due to the 
pollution of type II SNe, which explode on very short timescales (of the order of Myr to tenths of 
Myr) and produce roughly constant [cr/Fe] ratios. Then the [a/Fe] ratio decreases to the solar value 
in disk stars because of the bulk of Fe produced by type la SNe explode on timescales ranging 
from 30^0 Myr to 10 Gyr. This interpretation is known as "time-delay model" and represents the 
best interpretation of alpha-element abundance ratios in galaxies so far. Type II SNe originate from 
the core-collapse of massive stars (M > lOM©), whereas it is believed that type la SNe arise from 
mass-accretion in binary C-O white dwarfs systems, exploding when the white dwarf mass reaches 
the Chandrasekhar limit. Two main scenarios can lead to this situation: i) the single degenerate, 
where a white dwarf explodes after accreting material from a red giant or a main sequence com- 
panion (Whelan & Iben 1973); ii) the double degenerate, where the explosion follows the merger 
of two C-O white dwarfs (Iben & Tutukov 1984), due to loss of angular momentum from gravita- 
tional wave emission. There is not yet consensus on which scenario is better or if both are at work 
in galaxies. In both scenarios type la SNe start exploding after 30^0 Myr and continue to do so for 
a Hubble time. The maximum type la SN rate depends on the assumed progenitor model but also 
on the assumed star formation history, being short in systems with a strong burts of star formation 
(e.g. 0.5 Gyr in the galactic bulge), 1.0-1.5 Gyr in the solar vicinity and several Gyrs in slow star 
forming systems, such as irregular galaxies (Matteucci & Recchi, 2001). Therefore, the bulk of Fe 
is always produced with a delay relative to the a-elements produced by type II SNe. Manganese 
is an iron-peak element produced both in type II and la SNe. Nucleosynthesis calculations for 
massive stars by Arnett (1971) and Woosley & Weaver (1995; hereafter WW95), among others, 
show increased yields of Mn with increasing metallicity (actually with increasing neutron excess). 
Metallicity-dependent yields for type la SN have been largely unexplored; however, Ohkubo et al. 
(2006) showed that the yields of both Ni and Mn increase with metallicity in type la SNe, although 
this conclusion is sensitive to the details of the adopted hydrodynamical models. 



In this paper, we adopt three detailed chemical evolution models for the Galactic bulge, solar 
vicinity and Sgr dSph. In all the models we will adopt the same yields but different star formation 
histories tuned to reproduce the majority of the observational constraints on the three systems. We 
will show that to explain the [Mn/Fe] simultaneously in the three systems, the time-delay model 
alone is not enough, but we have to assume that the yields of Mn from type la SNe are metallicity- 
dependent. The paper is organized as follows: in section 2 we briefly describe the observational 
data, in section 3 we describe the chemical evolution models. In section 4 the nucleosynthesis 
prescriptions are given. In section 5 we discuss our results and in section 6 we draw our conclusions. 
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2. Observational data 

The observational data used in this paper for the Bulge and for the Sgr dSph are taken from 
McWilliam et al (2003). For Sgr dSph a part of the data come from the work of Bonifacio et 
al. (2000). Concerning the Bulge we added the mean values from three bulge globular clusters 
measured by Alves-Brito et al. (2005), Alves-Brito et al. (2006) and Ramirez & Cohen (2002). 
However, it should be noted that to compare the chemical abundances of globular clusters to the 
results of the chemical evolution model of the Bulge is not entirely safe; in fact, they might have 
had a different evolution when compared to the evolution of the Bulge or the Halo; moreover, it is 
still uncertain whether these globular clusters really belong to the Bulge system. 

The data for the solar neighbourhood shown in this paper come from the abundances reported 
by Reddy et al. (2003) and Reddy et al. (2006); they are not the only data available but we selected 
them in order to have a very homogeneous sample. 

3. The chemical evolution models 

For studying the chemical evolution of the solar neighbourhood we adopted the model of Frangois 
et al. (2004), which is an implemented version of the original model by Chiappini et al. (1997). This 
model assumes that the Galaxy formed by means of two main accretion episodes, one giving rise 
to the halo and thick disk and the other forming the thin disk. The infalling gas is always assumed 
to be of primordial composition. Detailed nucleosynthesis from low and intermediate mass stars, 
type la and type II SNe is taken into account. The IMF is taken from Scalo (1986). For details of 
this model see Frangois et al. (2004). 

For the Bulge we adopted the model of Ballero et al. (2007), which assumed a rapid formation 
timescale, of 0.3-0.5 Gyr, from gas accumulated during the Halo collapse. The efficiency of star 
formation (star formation per unit mass of gas) is 20 times higher (i.e. 2QGyr^^) than in the solar 
vicinity (IGyr^^). The IMF is flatter than in the solar vicinity, as required by the observed Bulge 
stellar metallicity distribution. See Ballero et al. (2007) for details of this model. 

To study Sgr dSph we have adopted the model developed for this galaxy by Lanfranchi et al. 
(2004; 2006). This model assumes that Sgr dSph formed stars over several Gyrs, as suggested 
by its color-magnitude diagram, and that a strong galactic wind, triggered by supernovae, was 
responsible for the gas loss. While the star formation efficiency is similar to the solar vicinity, the 
effect of significant gas loss is similar to a lower star formation efliciency. In particular, we assume 
that the star formation rate is given by a simple Schmidt law: 

l/r(f) = vo-gas (1) 

where ij/(t) is the star formation rate, v - IGyr"' is the star formation efficiency, namely the star 
formation rate per unit mass of gas, and o-gas is the surface gas density. The wind rate is then 
assumed to be: 

W(t) = Am (2) 

with A being a parameter whose best value for Sgr dSph is A = 13.0. The assumed IMF is also 
flatter than in the solar vicinity (Salpeter 1955). 

According to the time delay model, the predicted [a/Fe] ratios behave differently in the three 
systems because of the different star formation histories (see Matteucci 2003), in the sense that we 
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expect a long plateau with supersolar [a/Fe] ratios in the Bulge because of its fast star formation. 
In fact, in this case a lot of Fe is produced by type II SNe before the bulk of it is restored by 
SNe la. The opposite occurs for slow star formation systems, such as dwarf spheroidals. In this 
case the supersolar [a/Fe] plateau is much shorter than in the solar vicinity since when the type la 
SNe restore the bulk of Fe, the [Fe/H] in the interstellar medium is still quite low. This produces a 
situation where we expect low [a/Fe] ratios at low metallicities as opposed to the high [a/Fe] ratios 
at high metallicities expected for the Bulge. The observed behaviour of the [a/Fe] ratios in these 
three systems is in good agreement with these predictions (e.g. Francois et al. 2004; Ballero et al. 
2007; Lanfranchi et al. 2006). 

In the following sections we will test whether a time-delay scenario, where type la SNe produce 
the bulk of the Mn (as suggested Gratton 1989) can predict the observed [Mn/Fe] ratios in these 
systems. 

4. Nucleosynthesis Prescriptions 

As already mentioned, Fe and Mn are produced by type la and II SNe. In our chemical evolu- 
tion calculations we explore the consequences of three nucleosynthesis prescriptions: 1. fixed, 
metallicity-independent, Mn and Fe yields for type la and II SNe, 2. metallicity-dependent type 
II SNe yields with fixed, metal-independent, yields for type la SNe, and 3. metal-dependent Mn 
yields for both type la and type II SNe. 

In our first model we utilise the metal-independent yields from Frangois et al. (2004). For 
type II SNe they found that the WW95 results provide the best fit to solar vicinity data; in fact, 
no modifications were required for the yields of Fe computed for solar chemical composition. 
However, in order to fit the observed [Mn/Fe] trend Francois et al. modified the WW95 Mn yields 
for type II SNe, with a ~70% increase for 13-18 M© progenitors and a ~60% decrease for 30-40 
Mq progenitors. For type la SNe Frangois et al. (2004) employed the, unmodified, metallicity- 
independent Iwamoto et al. (1999) yields. Effectively, the Francois et al. (2004) used the age- 
metallicity relation and modified the Mn yields as a function of supernova progenitor mass to 
reproduce the observed trend of [Mn/Fe] with metallicity in the solar vicinity. Thus, these yields 
are fine-tuned to the formation timescale of the thin disk. We employed these yields for our initial 
chemical evolution calculations to predict the [Mn/Fe] trend for the metal-independent time-delay 
scenario. 

For our second nucleosynthesis prescription we chose to use the metallicity-dependent Mn and 
Fe yields of WW95 for the type II SNe, and we adopted the metallicity-independent theoretical 
yields by Iwamoto et al. (1999) for the yields from type la SNe. 

Finally, for the third nucleosynthesis recipe we again employed the WW95 metal-dependent 
yields for type II SNe but we also used metal-dependent Mn yields for type la SNe. For the metal- 
dependent type la SN Mn yields we modified the values of Iwamoto et al. (1999) by introducing 
the following metallicity dependence: 

/ >0.65 

Ym„(z)^Y'J:"""''\^-^J (3) 

where Z is the global metal content of the type la SN systems at birth. The form of the metal- 
dependent yields was obtained by requiring that the model predictions fit the observational data. 
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Very recently Badenes, Bravo & Hughes (2008) found that the Mn/Cr ratio in the ejecta of type la 
SNe depends on Z, and their modelled metallicity-dependentMn yield is similar to that in Equation 
(1), thus supporting our choice. Moreover, nucleosynthesis calculations by Ohkubo et al. (2006) 
suggested a dependence on metallicity of the yields of Ni and Mn in type la SNe. 
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Fig. 1. The results with Fran9ois et al (2004) set of yields. In Fig. the [Mn/Fe] vs [Fe/H] ratios in 
three different stellar systems; from left to right: Solar neighbourhood. Bulge and Sagittarius Dwarf 
Spheroidal. For the data see Sect. |2] 



5. Results 

5.1. Results with Frangois et al. yields 

In Figure 1 we show the results obtained with the yields of Frangois et al. (2004), and as expected, 
the data relative to the stars of the solar vicinity are well reproduced. Here and in the following 
figures, the results of all the models are normalized to the predicted solar abundances, namely 
the abundances predicted for the interstellar medium in the solar vicinity 4.5 Gyr ago. However, 
it should be noted that our predicted solar values for most of the considered abundances and, in 
particular for those studied in this paper, differ from the Asplund et al (2005) solar values by a 
maximum of ~ 0.06 dex. 

On the other hand. Figure 1 shows an unsatisfactory comparison between the time delay model 
predictions and the observed Bulge [Mn/Fe] trend with [Fe/H]: there is an ofi'set between the data 
and the results of the Bulge model, which predicts lower values of the [Mn/Fe] ratio for all metal- 
licities. 

Under the constraint of metallicity-independent Mn yields the time-delay model for the solar 
neighbourhood requires that the bulk of the Mn is produced by type la SNe. Because of the rapid 
Bulge formation (as suggested by by the alpha-element abundances), the contribution of type la 
SN to the bulge composition is smaller than the solar vicinity; therefore, the model predicts less 
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Mn production in the bulge in this paradigm, and thus lower [Mn/Fe] ratios (contrary to the obser- 
vations). 

Moreover, with this set of yields the results of Sgr dSph model also do not match the observed 
[Mn/Fe] trend with [Fe/H]: the theoretical curve passes through the observed [Mn/Fe] values only 
at [Fe/H]~ -1.5 dex, but for higher metallicities the predicted trend is too high compared to the 
data. This is understood as a result of the high frequency of type la SNe in the the Sgr dSph model 
(suggested by the low alpha/Fe ratios), which produce more Mn for metallicity-independent Mn 
yields in the time-delay scenario. Thus, the prediction of this paradigm is for enhanced [Mn/Fe] 
ratios in the Sgr dSph (again contrary to the observations). 
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Fig. 2. The results from the WW95 metallicity-dependent type II SN yields, with constant, Iwamoto 
et al. (1999) type la SN Mn yields. The [Mn/Fe] versus [Fe/H] ratios in three different stellar 
systems, from left to right: Solar neighbourhood. Bulge and Sagittarius Dwarf Spheroidal galaxy. 
For the data see Sect. |2] 



5.2. Results with Woosley & Weaver yields 

To improve the match between our models and the data, in particular for the Bulge and Sgr dSph, 
we now explore the use of the metallicity-dependent WW95 yields from type II SNe, but with 
the constant (i.e. metallicity-independent) Mn and Fe yields for type la SNe from Iwamoto et al. 
(1999); the results are shown in Figure 2. 

The results of the solar neighbourhood model still follow nicely the observational data. 

The model results for the Bulge are also more promising; in fact, up to solar [Fe/H] the pre- 
dictions well trace the observational data. However, for supersolar metalhcities the results of the 
model tend to shghtly underestimate the value on [Mn/Fe] in the Bulge stars. Because of the high 
star formation rate in the Bulge type la SNe occur at late times (and high [Fe/H]). Thus, the Mn 
deficit at supersolar metallicities in this Bulge model might reasonably be due to a problem with 
Mn from type la SNe. 
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The results from this model for the Sgr dSph are still in poor agreement with the data. The 

model shows a strong increase of the [Mn/Fe] ratio at [Fe/H] 1 .5 dex, at variance with the data. 

Given the prominent role of type la SNe in Sgr dSph this Mn over-prediction is likely related to the 
type la SN Mn yields. 
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Fig. 3. The results with metallicity dependence in both type II SNe and type la SN yields for Mn. 
In Fig. the [Mn/Fe] vs [Fe/H] ratios in three different stellar systems; from left to right: Solar 
neighbourhood, Bulge and Sagittarius Dwarf Spheroidal. For the data see Sect.|2] 



5.3. Results with Woosley & Weaver yields and metallicity dependent Mn yields from type la 
SNe 



In this last set of chemical evolution models both the type la and II SN yields of Mn are dependent 
on the metallicity. In particular, the Mn yields from type la SNe follow Equation (1). In this case the 
yields are constrained to maintain agreement between the chemical evolution model predictions for 
[Mn/Fe] versus [Fe/H] and the solar vicinity observational data. At the same time, the supersolar 
data in the Bulge are nicely fit. Thus, the use of this new set of yields had improved the predictions 
of the Bulge model, presumably because the Mn contribution from late time, high metallicity, type 
la SNe are more coiTectly modelled. Finally, the most interesting result is the excellent agreement 
of the Sgr dSph model with the data. This is due to the reduced contributions of the Mn enrichment 
by low metallicity type la SNe, compared to the time-delay metallicity-independent model. 

A very important conclusion is that the yield of Mn produced by type la SN must be metallicity- 
dependent in order to fit the Sgr dSph abundance data. Moreover, we get excellent agreement also 
for the Bulge and solar vicinity stars. It is worth emphasising that it was not possible to obtain 
this result without the simultaneous comparison of the three independent stellar systems, charca- 
terised by different star formation timescales. This hypothesis, that the type la SN Mn yield must 
be dependent on metallicity, was first suggested by Mc William et al. (2003). 
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6. Conclusions 

We show quantitatively, for the first time, that galactic chemical evolution models demand that the 
type la SN Mn yields must be metallicity-dependent, by simultaneously fitting the [Mn/Fe] trend 
in three independent stellar systems: the Galactic bulge, the solar vicinity and the Sgr dSph. We 
have adopted detailed chemical evolution models, chosen to reproduce the main chemical features 
of these three stellar systems, using same element yields but different star formation histories. 

We have shown clearly that the "time-delay model" alone cannot explain the different behaviour 
of the [Mn/Fe] in Bulge, solar vicinity and Sgr dSph, as it does for [c/Fe] ratios. In fact, in the 
case of the [a/Fe] ratios, the same models can well reproduce the observed [a/Fe] versus [Fe/H] 
trends in the three stellar systems without invoking metallicity-dependent yields. This is because 
the production of a and Fe elements is only mildly dependent on the initial stellar metallicity. We 
propose a functional form of the metallicity-dependent Mn yields in type la SNe (oc Z"^^) which 
should be tested by nucleosynthesis calculations. Our result is supported by Ohkubo et al. (2006), 
who suggested that the yields of Mn from type la SNe have a strong dependence on metallicity, 
and by a recent finding by Badenes et al. (2008), who showed that the Mn to Cr ratio in the ejecta 
of type la SNe depend on the metallicity of the progenitor star in a very similar way to what we 
have suggested. 

It is worth noting that in this work we have taken the Reddy et al. (2003, 2006) abundances for 
the solar vicinity stars at face value. However, we note that McWilliam et al. (2003) found that a 
zero-point correction of +0.10 dex to the Reddy et al. (2003) [Mn/Fe] values was required if the 
sun is typical of stars in the range -0. 10 < [Fe/H] <i +0. 10. We have re-investigated this issue and 
found that together the Reddy et al. (2003) and Reddy et al. (2006) data indicate a mean [Mn/Fe] 
of solar [Fe/H] stars at -0. 10+0.01 dex, with an rms scatter about the mean of 0.036 dex. Since the 
solar Mn/Fe ratio is 2.77 sigma above the mean for solar [Fe/H] stars (equal to 99.7% above the 
mean in a gaussian distribution), this zero-point shift is most likely due to systematic error in the 
Reddy et al. (2003, 2006) results, rather than a real Mn enhancement in the Sun. This is supported 
by the fact that other abundance studies require quite different zero-point corrections; for example 
the differential study, relative to the sun, by Feltzing & Gustafsson (1998) required a zero-point 
[Mn/Fe] correction of -0.06 dex. 

Unfortunately, we are unable to make as good a fit to the measured [Mn/Fe] trends if we employ 
a 0.10 dex zero-point shift to the Reddy et al. (2003, 2006) data. To do so would require a more 
complicated form of the metallicity-dependent Mn yields from type la SNe, instead of the simple 
power-law form adopted here. Alternatively, other complications might also be responsible, such 
as the type II SNe yields of WW95, or variations in the star formation and/or the mass-loss rates 
for Sgr dSph. Rather than explore a myriad of possibilities that might give a perfect fit to all the Mn 
trends, we prefer to focus on our main quantitative result: that the trend of [Mn/Fe] with [Fe/H] in 
the Galactic disk, the Bulge, and the Sgr dSph can only be understood if the yield of Mn from type 
la SNe is metallicity-dependent. 

Finally, we note that the low [Mn/Fe] values seen in Sgr dSph stars near solar metallicity is due, 
at least in part, to the lower metallicity of old Sgr dSph stars compared to the solar neighborhood. 
It is the low average [Fe/H] of the type la SNe in the Sgr dSph that results in lower [Mn/Fe] values, 
thanks to the metallicity-dependence of the Mn yield. Thus, the [Mn/Fe] ratios of solar metallicity 
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stars in dwarf galaxies is linked to the mean metallicity of these galaxies. Since the metallicity 
distribution function of stars in dwarf galaxies is sensitive to the mass outflow and star formation 
rates, the value of these parameters should be constrained by the observed [Mn/Fe] ratios. Indeed, 
we would expect that the more massive dwarf galaxies (e.g. LMC) would show a similar effect 
on [Mn/Fe], but less extreme, than in lower mass dwarf galaxies; thus, there should be families of 
low-mass galaxies with different [Mn/Fe] trends. 
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